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The semileptonic decays of fie and Q.b are treated in the framework of a constituent quark model 
developed in a previous paper on the semileptonic decays of heavy A baryons. Analytic results 
for the form factors for the decays to ground states and a number of excited states are evaluated. 
For fib to fic the form factors obtained are shown to satisfy the relations predicted at leading 
order in the heavy-quark effective theory at the non-recoil point. A modified fit of nonrelativistic 
and semirelativistic Hamiltonians generates configuration-mixed baryon wave functions from the 
known masses and the measured Aj — > Ne'^v rate, with wave functions expanded in both harmonic 
oscillator and Sturmian bases. Decay rates of fii, to pairs of ground and excited Q.c states related 
by heavy-quark symmetry calculated using these configuration-mixed wave functions are in the 
ratios expected from heavy-quark effective theory, to a good approximation. Our predictions for the 
semileptonic elastic branching fraction of fig vary minimally within the models we use. We obtain 
an average value of (84± 2%) for the fraction of 0,^ — > H'-*' decays to ground states, and 91% for 
the fraction of Q,c fi'*' decays to the ground state fi. The elastic fraction of fif, ^ Q.^ ranges 
from about 50% calculated with the two harmonic-oscillator models, to about 67% calculated with 
the two Sturmian models. 
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I. INTRODUCTION AND MOTIVATION 



Semileptonic decay of hadrons are of interest for two basic reasons; they are the primary source of information 
for the extraction of the Cabbibo-Kobayashi-Maskawa (CKM) matrix elements of the Standard Model from 
experiment, and the study of the semileptonic decays of baryons provides information about their structure. In 
this manuscript, we present results of a calculation of the form factors and rates of the semileptonic decays of 
heavy 11 baryons (fig) obtained using a constituent quark model. This work is similar to a recently published 
calculation of the semileptonic decays of heavy A baryons (Ag) (hereinafter referred to as I). Although our 
motivation for the present work is similar, it is briefly recapitulated here for completeness. 

The heavy quark effective theory (HQET) 3| has been a very powerful tool for the extraction of the CKM 
matrix elements from data on the semileptonic decay of mesons, especially for the decays of heavy mesons to 
heavy mesons where a number of relations that simplify this extraction are provided. For the semileptonic 
decays of a heavy baryon to another heavy baryon, HQET makes predictions that are analogous to those made 
for heavy-to-heavy meson decays: the six form factors that describe the decays to the ground-state heavy 
baryons are replaced by fewer form factors called Isgur-Wise functions; the normalization of at least one of 
these Isgur-Wise functions is known at the non-recoil point; corrections to this normalization first appear at 
order l/mq; and corrections can be systematically estimated in a I/itiq expansion. Note that for baryons, 
the number of Isgur-Wise functions needed at leading order depends on the flavor-spin structure of the parent 
and daughter baryons. For the semileptonic decays Qq — > Qq, only two functions are required to describe the 
semileptonic decay in the heavy-quark limit. In the case of a heavy baryon decaying to a light baryon, HQET 
makes predictions that are not as powerful as in the heavy-to-heavy case. For example, for the semileptonic 
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decays Qq — > fl, the leading-order HQET prediction is that the number of independent form factors decreases 
from six to four for a daughter ft with spin 1/2. 

While HQET has been tremendously successful and useful in treating the semileptonic decays of heavy 
hadrons, it has limitations. It is a limit of QCD that applies only to hadrons containing heavy quarks, and 
for the decays of such hadrons, it only predicts the relationships among form factors, not their kinematic 
dependence. In addition, the predictions of HQET are valid only as long as the energy of the daughter hadron 
is much smaller than the mass of the heavy quark. These limitations mean that the predictions of HQET must 
be augmented by information arising from other approaches to hadron structure. 

While some work has been done in modeling the form factors for the semileptonic decays of heavy baryons, 
to the best of our knowledge little has been done in treating the decays to excited baryons. In I, the existing 
models of semileptonic decay in both meson and baryon sectors have been discussed. We have also outlined 
the procedure used to obtain the form factors and decay rates for semileptonic decay of Ag baryons using a 
non-relativistic quark model. Here we will focus only on the new aspects relevant to the semileptonic decay of 

Very little theoretical or experimental work has been published to date on the semileptonic decay of fig 
baryons. Boyd and Brahm Q have used HQET to show that the fourteen form factors that describe the decays 
fih — > ^l^eD and fif, — > 0*(J^ — 3/2^)eD can be parametrized in terms of two nonperturbative functions at 
leading order, and in terms of five additional nonperturbative functions and one dimensional constant at order 
l/iric- A l/rric expansion of the form factors for fif, fic*^ has also been carried out by Sutherland 0|, who 
treated the effect of the fig — ilg mass splitting on the form factors evaluated at first order in l/rric- A Bjorken 
sum rule for the semileptonic decays of fh, to the ground state and low-lying negative-parity excited-state 
charmed baryons has been derived by Xu [3, again in the heavy quark limit. To the best of our knowledge, 
there are no calculations performed outside of the HQET framework, which motivates the present study. 

There has recently been some progress in experiments dealing with these decays. The CLEO-c collaboration 
has published evidence for the observation of the decay il^e+i^, and have measured the product of the 

branching fraction and cross section to be i?(ri° Vl~e^v) ■ a{e^e^ ^cX) = 42.2 ± 14.1 ± 11.9 fb. The 
ARGUS JJ and BELLE ^] collaborations have also seen evidence for the Q,^e^v decay, but no quantitative 

value for the branching fraction has yet been published. 

The procedure we follow here to calculate the form factors and rates for VIq semileptonic decay is similar to 
that used for Ag decays in I. In a quark model context the flavor part of the Ag wave function is anti-symmetric 
under exchange of quarks 1 and 2, which comprise the light diquark system. This requires the spin-momentum 
part of the Ag wave function to be anti-symmetric under exchange of the two light quarks to maintain the 
appropriate Pauli exchange symmetry. On the other hand, fig baryons are described by a symmetric flavor 
wave function for the light diquark system. As a result, the spin-momentum part of the wave function has to 
be symmetric. This basic difference between the wave functions of Ag and f2 makes calculations of the form 
factors and the corresponding decay rates for their semileptonic decays significantly different. 

This manuscript is organized as follows: in Section II we discuss the hadronic matrix elements and decay 
rates. Section HI presents a brief outline of heavy quark effective theory as it relates to the ilg decays that 
we discuss. In Section IV we describe the model we use to obtain the form factors, including some description 
of the Hamiltonian used to generate baryon wave functions. Our analytic results are discussed and compared 
to HQET results in Section V, our numerical results are given in Section VI, and Section VII presents our 
conclusions and an outlook. A number of details of the calculation, including the explicit expressions for the 
form factors, are shown in three Appendices. 

II. MATRIX ELEMENTS AND DECAY RATES 

The transition matrix element for the semileptonic decay f2g — > ^qH-vi^ is written 

T^^VQ,u-a^{i--f,)u,,{ii,{p',s')\j^\nQ{p,s)), (1) 

where Gp/^f^ = 5^/(8-/^/^) is the Fermi coupling constant, Mw is the intermediate vector boson mass, Vqq 
is the CKM matrix element, 1(77^(1 — 75)^1/, is the lepton current, and = 97^^(1 — 75)Q is the left handed 
current between quarks Q and q. The hadronic matrix element of J^^ is described in terms of a number of form 
factors. 
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For transitions between ground state (J^ = 1/2+) baryons, the hadronic matrix elements of the vector 
(Vfj, = q^fxQ) and axial (A^ = 97^75(9) currents are 

{B,{p\s')\VMp,^))=W,s') (F,{q^)^^ + F,{q^)^+Fs{q^)^)u(p,s), (2) 

V mBQ rUB, ) 

(B,(p',s')l^jSQb,s)) ^u{v'.s') ('Gi(q2)^^ + G2(g2)i + G3(g2)^)75^^(p,s), (3) 

V mBQ / 

where the Fi and Gi are form factors which depend on the square of the momentum transfer q = p — p' between 
the initial and the final baryons. Similar expressions for transitions between = 1/2"*" ground states and final 
state baryons with other spins and parities are given in I. If J > 3/2 these expressions involve a fourth pair of 
form factors F4 and G4. 

Expressions for the differential decay rates for semileptonic decays both including and ignoring the mass of 
the final leptons are given in I. These expressions can be integrated to yield the decay rates reported later in 
the present paper. 

III. HEAVY QUARK EFFECTIVE THEORY 

In most applications of HQET, the aim has been to constrain the hadronic uncertainties in the extraction of 
CKM matrix elements such as Vub and Vcb- In this section, we take a different tack; we examine the predictions 
of HQET for decays of a heavy f2 into a number of the allowed excited heavy daughter baryons, with the aim 
of comparing these predictions with the form factors that we obtain in our model. 

A. Structure of States and Parity Considerations 

In a heavy excited baryon, the light quark system has some total angular momentum j, so that the total 
angular momentum of the baryon can be J = j ± 1/2. These two states are degenerate because of the heavy 
quark spin symmetry. It is useful to show explicitly the representation we use for these two degenerate baryons. 
In the notation of Falk |9|, we write u'^-Xi^ = u^'^-'^' W) - W)^ with 

^t^i-H ^ {v')uq{v'). (4) 

Here, uq(v) is the spinor of the heavy quark, with v being the four- velocity and A^^^---^^ {v') is a tensor that 
describes the spin-j light quark system. This tensor is symmetric in all of its Lorentz indices, meaning that 
u'^i - ''^ (ti') is also symmetric in all its Lorentz indices. Both u^j^-i^^ {v') satisfy the conditions 
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0, .g^,^,«'^-''^(i;) = 0, (5) 



where and /i/ indicate any pair of the indices /ii . . . /ij. The state with J = j + 1/2 also satisfies 

= 0. (6) 

Further details of the structure and properties of these tensors are given in Falk's article . 

At this point, it is useful to discuss the parity of the states, which is determined by the parity of the light 
component. A spin-j light quark component with parity (—1)-' is said to have 'natural' parity, unnatural parity 
otherwise. The natural-parity light quark systems therefore have = (271)"*" or = (2n -I- 1)~, with n a 
positive integer or zero. The natural-parity light quark systems are represented by tensors, while those with 
unnatural parity are represented by pseudo-tensors. Since the parity of the baryon is that of the light quark 
system, we may refer to the baryons as being tensors or pseudo-tensors, with the understanding that this really 
refers to the light-quark component of the baryon. It is thus convenient to divide the decays we discuss into 
two classes, those in which the daughter baryons are tensors, and those in which they are pseudo-tensors. 
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B. Heavy to Heavy Transitions 



First, we note that the ground state of the Qq has a symmetric flavor wave function for the hght diquark, and 
so has a spin wave function that is also symmetric, corresponding to a total spin equal to one in the light quark 
component of the wave function. This state is therefore the spin-1/2 member of the lowest-lying (1/2+, 3/2+) 
multiplet. The Falk representation of this state is 



(7) 



where u{v) is a Dirac spinor. This state is a pseudotensor, and we begin with a discussion of decays to other 
pseudotensor states. We are interested in the matrix element 



A=<n[*\v\j)\cTb\nb{v) >, 



(8) 



where c and b are the heavy quark fields, and F is an arbitrary combination of Dirac matrices. With the use of 
HQET, we may write this matrix element as 



(9) 



to leading order. Here, M^^^^^^.^, is the most general tensor that we can construct, fi" is the Falk representation 
of Eq. iQ, and u'^^'"'^^ {v') is the analogous representation of the daughter baryon. M^^,,,!^^-^ may not contain 
any factors of w^. or gfj^fij , and therefore takes the form 



M, 



(10) 



Thus, two independent form factors, ri[''l{v ■ v') are needed to this order, regardless of the spin of the final 
baryon. 

Applying these results to the specific case of = 0^, we find, for = 1/2^ 



w - 



(0), 



^/3 



w + 1 (0) 



V3 



(11) 



where w — v ■ v' . In this case, the daughter baryon is a singlet, and has no Lorentz indices. This means that 
the term in g^i, is not present, and only the term in 772 contributes to the matrix element. 
When jP = 1+, we find, for = 1/2+, 



and for 

Gi = 



Fi 
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= -G4 = 


1 

"71 



















(12) 



'V3 



m 



(13) 



In these two sets of equations r][^l are two universal functions of the Isgur-Wise type. We note that there exist 
several multiplets that have the same quantum numbers as this ground state multiplet. The same is true in 
the case of Aq baryons, and these excited Aqs were identified as radial excitations of the ground state. In the 
case of the fig, such baryons are indeed excitations of the ground state, but they are not necessarily radial 
excitations. Some of these excitations are orbital excitations. However, independent of whether the daughter 
baryon belongs to the ground state (1/2+, 3/2+) or one of the excited multiplets, the expressions above for the 
form factors are valid. The explicit forms of the Tyl^-* will depend on the details of the structure of the daughter 
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baryon. For the ground state we know that at the non-recoil point, r][^\v.v' = 1) = —1, while the normalization 
of 7^2^'' is not known. The negative sign of the normalization of r^j^^ arises because we have chosen a positive 
sign for the g^^„ term in Eq. ((Tn|l . 

For = we find for = 3/2-, 



F3 
Gi 
G3 



(2w - 1)77^^^ +2{w^ - 1) 77^- 

(2) _L0^„_ 1^.(2) 



(2) 



15 



(2) 



(2) 



1 



30 
2 



(2w + l)77f ^ + 2 (w^ - 1) r]^ 



(2) 



, G2 = -2\/ — 



77f^ + (z« + l)r;f 



77^^ + 2 (zi; + 1) ryf 1 , G4 = + 1)77}- 

J V 15 



(2) 



and for = 5/2 

Fi = 
Gi = 



1 
1 

7! 



(2) I / 1 ^ I- 
rjl' + {w-l) ri'^ 



(2) 



F2 — G2 — 0, -Fa 



-G3 



V3 



ri(2) p 



(2) 



77p^ + (ii; + 1) 7;^- 



(2) 



As with the previous example, the functions rjl ^ ^re Isgur-Wise form factors common to both decays. 
For f = 3+, wc find for = 5/2+, 



Fi = 

Fz = 

Gi = 

G3 = 



1 



3^7 
2 



(3w - 2)7;f ' + 3 (tw^ - 1) 77^' 



(3) 



F2 = -2^/- 



77^^ + (tw - 1)77^' 



(3) 



3\/7 
1 



77f)+3(u;-l)7yf ] , F4 = -^(77; - 1)77^ \ 



3\/7 
2 

3\/7 



(37« + 2)7;f' + 3 (77;2 - 1) 



(3) 



G2 = -2j3 



77f^ + (7i; + 1)77^^ 



77f ^ + 3 (7« + 1) 77^^^ 



G4 — 



3\/7 



{w + l)77i' 



(3) 



and for = 7/2+ 
Fx = 
Gi = 



1 
1 

71 



(3) I / IN (3) 

77i ^ + - 1) 77^ ^ 

77P' + (W + 1) 77^^ 



, -F2 — G2 — 0, F3 



-G3 



V3 



F4 



V3 



(3) 



(14) 



(15) 



(16) 



(17) 



The functions ri\ ^ are Isgur-Wise form factors common to both decays. The normalizations oi r\\2 are not 
known. 

For the tensor decays, the matrix element again takes the form shown in Eq. ((SJ, but M^j...^^,y must now be 
a pseudotensor. The only form that we can write is 



fj.i...fj.ji' 



(18) 



When applied to the 1/2+ singlet daughter baryon, there is no way to create this pseudotensor, so such am- 
plitudes vanish at leading order. For the other spin states, after some manipulation, we can express the form 
factors in terms of the set of Isgur-Wisc functions t'^^\w). 
For = f-, wc find for the 1/2^ state 



-Fi — Gi — 0, F2 — F3 — — G2 — G3 



(19) 
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while for 3/2 state, the form factors are 

F2 = G2 = 0, G3 = -F3 = 2Fi = -2Gi = -^r(i\ F4 = -^{w- 1)t^^\ G4 = ^{w + 1)t^^K (20) 

v3 v3 v3 

For = 2+, starting with 3/2+, the form factors are 



Gi = ^(l + «,)r(2), G3 = -^/jH + 2)r(^^- (21) 
V oU V 10 



For the 5/2+ state, the form factors are 

F2 = G2 = 0, Gs = -Fs = 2Fi = -2Gi = ^^^'^^"^^ " 7!^"" + ^^^^ 

The normalizations of none of the r'*' are known. 

We do not present the predictions for the decays of fig to light O states, as the HQET predictions are not 
as useful as they are in the case of heavy to light Ag decays. For instance, the decays of the ground state fig 
to an f2 with = 1/2+ are described in terms of four form factors in HQET, instead of six in general. While 
this small simplification is no doubt useful we do not pursue it here. 

IV. THE MODEL 
A. Wave EYinction Components 

In our model, a baryon state has the form 

\Aq{p,s)) = 33/* y rf=^PprfVC^^iQkl(Pl,Sl)92(P2,S2)93(P3,S3)), 

where "^(Pi ^ P2) = -^(Pi +P2 — 2p3) are the Jacobi momenta, is the totally antisymmetric 

color wave function, and "^^q ~ ^Aq^^AqXAq is a symmetric combination of flavor, momentum and spin wave 
functions. The flavor wave functions of fig and S are 

(jjfiQ — ssQ, (f>so = ssu, (ps- = ssd, 

which are symmetric in quarks 1 and 2. The momentum-spin parts of the wave functions must therefore be 
symmetric in quarks 1 and 2 to keep the overall symmetry. The symmetric spin wave function xf/2> ^^'^ 
mixed symmetric spin wave functions Xi/2' X1/2 usual eigenstates of total spin made of three spin-1/2 

quarks. 

The momentum wave function for total L = £p + £x is constructed from a Clebsch-Gordan sum of the wave 
functions of the two Jacobi momenta Pp and p^, and takes the form 

■^LMn.e^nxexiPp^Px) = X^(-^-^Kp"Z, ^A^f - TO)V'np<!pm(Pp)V'n^^M-m(PA)- 

m 

The momentum and spin wave functions are then coupled to give symmetric wave functions corresponding to 
total spin J and parity (— 1)(V+'a)j 

*JM = J2^JM\LMl,SM - ML)tpLMi^n,e,n>,edPp^P\)xs{M - Ml) 

Ml 

= [tpLMLn,e,n^ex{Pp,Px)Xs{M-ML)]j^j^. (23) 
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The full wave function for a state A is built from a linear superposition of such components as 



"fAJPM ^(bAY.rif'i'jM. (24) 

i 

Here (f>A is the flavor wave function of the state A, and the rjf are coefficients that are determined by diagonalizing 
a Hamiltonian in the basis of the '^^j^^- For this calculation, we limit the expansion in the last equation to 
components that satisfy N < 2, where N = 2 (rip + n\) + £p + £\. Consistent with this is the fact that the states 
we discuss all correspond to < 2. 

The wave functions for fig with = 1/2+ have the form 



T i 

*l/2 + A/ = Wc; 



??! '^■0000000 (Pp, Pa) + % '''0ooiooo(Pp,Pa) + % V'ooooio(Pp, Pa) Xi/2(^'^) 



+ f]r'''0OOOlOl(Pp,PA)x^/2(^^) IV'lAfi010l(Pp,PA)x?/2(*^- A^i) 



1/2, M 



1/2, Af 



+ V7 



'02Afi,O2Oo(Pp, Pa)X3/2(A^ " ^l) 

V'2Af^0002(Pp, PA)xf/2(^^ - Ml)] ^ J . (25) 

The complete expressions for fig wave functions of different spins and parities are given in Appendix ^ A 
simplified version of the model would truncate the expansion of the wave functions, giving 

*OQa/2+A./ = 0nV'OOOOOo(Pp,PA)XV2(^-^) (26) 

for the ground state. There are a number of daughter baryons that have an overlap with the ground state fig 
(in the spectator approximation that we use), even when we limit the discussion to states with N < 2. There 
are three states with — 1/2+, two with — l/2~, two with = 3/2~, four with = 3/2+, one with 
= 5/2~, two with = 5/2+, and one with — 7/2+, all of which occur in the N < 2 bands. The 
single-component representations of these states are 

*aQ,i/2+Af = 0nV'oooooo(Pp,PA)xv2(^)' 
*aQ,l/2+Af = '/'oV'ooooio(Pp,Pa)xv2(^^): 

*nQ,l/2+Af = </'o[V'2A/i0002(Pp,PA)xf/2(^^-^i)]l/2,Af, 
*OQ,l/2-Af = 0o['01A/iOOOl(Pp,PA)Xl/2(^^ - ^^i)]l/2,A/, 

*OQ,i/2-Af = (t>n[i'iMi^aoaiiPp,P\)xi/2iM ~ Ml)]i/2M, 

*Oq,3/2-A/ = 0o['01A/i,OOOi(Pp,Pa)Xi/2(^^ - ^^i)]3/2,Af, 

*OQ,3/2-Af = </'o[V'lA/i000l(Pp,PA)xf/2(^^-*^i)]3/2,A//, 

*OQ,5/2-Af = 0o['01Afi,OOOl(Pp,PA)xf/2(^^ - *^i)]5/2,Af, 

*aQ,3/2+Af = 0oV'OOOOOo(Pp,PA)xf/2(*^): 

*aQ,3/2+Af = 0oV'OOOOlo(Pp,PA)xf/2(^). 

*nQ,3/2+Af = '/'o[?A2A/x,0002(Pp,PA)xf/2(*^-^i)]3/2,M, 

*nQ,3/2+Af = </'o[V'2A/i0002(Pp,PA)XV2(^^-^^i)]3/2,A/, 

*nQ,5/2+Af = '/'o[?A2A/i,0002(Pp,Pa)XV2(^^ - ^i)]5/2,Af, 

*nQ,5/2+Af = '/'o[?A2A/x,0002(Pp,PA)xf/2(^'^ " AfL)]5/2,M, 

*nQ,7/2+Af = 0o[V'2A/iOOO2(Pp,PA)xf/2(^^ - ^^i)]7/2,A./- (27) 

A common choice for constructing baryon wave function is the harmonic oscillator basis. One advantage of 
using this basis is that it facilitates calculation of the required matrix elements. However, it leads to form factors 
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that fall off too rapidly at large values of momentum transfer. We therefore also use the so-called Sturmian 
basis 10] . In this basis, form factors have multipole dependence on g^, which is what is expected experimentally. 
The explicit wave functions in momentum space are 



2n\ 



1 



(28) 



in the harmonic oscillator basis, and 



{n + L + \)\ 



(f^ + l) 



L+2 



Pr 



2-^+2 J fp^-f3 



+ f3^ 



J^Lm(p) 



(29) 



in the Sturmian basis. The L'^{x) are generalized Laguerre polynomials and the Pn^'^^ (y) are Jacobi polynomials, 
with p = |p|. 



B. Hamiltonian 



The phenomenological Hamiltonian we use has the form 

3 3 



^ E + C.„ + i: ( - + S. ^ S,*3(,,) ) , ,30) 

1=1 i<cj—l 



with rij = \ri ~Vj\. The spin independent confining potential consists of a linear and a Coulomb component, 
and the spin-dependent part of the potential takes the form of a contact hyperfine interaction. Spin-orbit and 
tensor interactions are neglected. We note here that acoui, Q^hyp, b, Cqqq, and rrii are not fundamental, but are 
phenomenological parameters obtained from a fit to the spectrum of baryon states. In the sum Yl^i=i of the 
kinetic energies of the quarks, each term has either the usual non-relativistic form given by 

9 

2m,j 

or a semi-relativistic form given by 



K,^[m, + ^] , (31) 



K^^Jpi + mf. (32) 



C. Obtaining the Form Factors 



1. Q,Q Q,q 



Here, we illustrate the procedure we follow to obtain the form factors, using the decay of the Bq to the ground 
state Bq as an example. We note here that Bq represents JIq and Bq refers to any of the ^Iq or S in their ground 
state. We show the procedure only for the vector current matrix element from Eq. |(5J, with the assumption 
that the parent Bq is at rest and the daughter Bq has three momentum p. The left-hand side of Eq. (jSJ is 
evaluated using the quark model, after the operator = q^y^Q has been reduced to its Pauli (non-relativistic) 
form. Specific values for the index p, are chosen, as well as specific values of s and s' . By making three sets 
of such choices, three equations for the Fi in terms of the quark-model matrix elements of three operators are 
obtained. This system of equations is then solved to obtain the expressions for the form factors. In the specific 
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case at hand, choosing s = s' ~ +1/2 and n — 0, for instance, leads to 

{Bg{p,+)\q^oQ\BQ{0,+)) = I d'pyp\d'ppd^pxC^*C^n''^{+) 

X {qW2<l\qhoQ\qiq2Q)'i'%Qi+)- 

(33) 

where 

{qW2q\qhoQ\qiq2Q) - {qW2\qiq2){q\qhoQ\Q). (34) 

The matrix element ( 9192 19192) gives (5-functions in spin, momentum and flavor in the spectator approximation. 
Using the J-functions in momentum and flavor, the integral simplifies to 

(B,(p,+)|Oo|Sq(0,+)) = y'd3p,dVV'B,(p;„PA)-4^+s,(^?o)V's«(Pp,PA), (35) 



with p'p — Pp, p'^ = Pa — 2-\/3/2mcrp/TOB^, where nia- is the mass of the light quark and Oq = q^joQ- It is 
useful for us to define 



(Of.) = xiy)Ss',s,6,',sAqiP3,4)\0^\Q{P3,S3))xBc,is). (36) 



where P3 = P — y |pa and p3 — — y |pa- We have 

-4b^b,((^m) - ^(9(P3, t) \0^\Q{P3, T)) + ^(<Z(P3, i) \0^\ Q(P3, i)), (37) 

where both parent and daughter baryons are in the ground state. 

After the spin matrix elements are evaluated, the momentum integrals are performed using both bases for 
the momentum wave functions shown earlier. The analytic results for the form factors for the Qq — > fi^ decays, 
evaluated using the truncated basis of Eq. (|27|l . are given in Appendix^ For decays to excited states, the 
calculation of the form factors is a little more involved, but the basic idea is as outlined here. 



2. fie ^ n'*' 

The calculation of form factors for flc O*^*^ decays is similar to that described above, up to a question of 
symmetry. The flavor wave functions of flc and are — ssc and = sss. In a semileptonic decay process 
the charm quark of the Qc decays into an s quark, which can be any of the three s quarks of the il. Thus, we 
need to evaluate 

i= (^(f](p, +)|o^|f^Q(o, +)> + (o(p, +)\{u}OMo, +)) + (f^(p, +)\moMo. +))) ■ 

The factor comes from the normalization. The wave functions for the O states are fully symmetric un- 
der interchange of any of the quarks, so each of the permuted matrix elements reproduces the one without 
permutations. The result is that we must calculate 

V3{n{p,+)\o^\nQ{o,+)), 

to obtain the form factors for flc fl. 

The procedure described in this subsection is relatively straightforward to implement in the harmonic oscillator 
basis, largely due to the fact that the Moshinsky rotations have been treated by a number of authors, and are 
also fairly simple to calculate. In particular, the fact that the 'permuted' wave function can be written in terms 
of a finite set of transformed wave function components is another feature that makes the harmonic oscillator 
basis attractive for calculations like these. In the Sturmian basis, however, the permutation of particles requires 
an infinite sum of transformed wave functions. This sum could be truncated at some point in a calculation such 
as this. However, at this point we do not examine decays to daughter f2's in the Sturmian basis. 
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V. ANALYTIC RESULTS AND COMPARISON WITH HQET 

The analytic expressions that we obtain for the form factors are shown in Appendix^ for both the Sturmian 
and harmonic osciUator bases. The results shown there are valid when the wave function for a particular state 
is written as a single component, in either expansion basis. As mentioned earlier, one of the advantages of the 
Sturmian basis is that it leads to form factors that behave like multipoles in the kinematic variable, and this is 
seen in the forms that we display. 

At this point, it is instructive to compare, as far as possible, these analytic forms with the predictions of 
HQET. While HQET does not give the explicit forms of the form factors, a number of relationships among the 
form factors are expected, and any model should reproduce these relationships. In what follows, we restrict 
our comparison to the predictions that are valid at the non-recoil point, as we have ignored any kinematic 
dependence beyond the Gaussian or multipole factors shown in Appendix ^ In addition, we focus on the 
predictions for heavy to heavy transitions. 

The quark model states we use are constructed in the coupling scheme 

\j'',L,S>^\[{lph)L{si2S3)s]j), (38) 

where the notation {ab)c means angular momentum c is formed by vector addition from angular momenta a 
and b. The parity P is (—1)^''+'^^, the total spin of the two light quarks in the baryon is si2, and S3 is the spin 
of the third quark, taken to be the heavy quark. 

The HQET states are assumed to have the coupling scheme 

|J^,J>= |{[(VA)iSi2]^.S3}^), (39) 

where j is the total spin of the light component of the baryon, so that J = j zL 1/2. The states of one coupling 
scheme are linear combinations of the states of the second. In particular, we find 



1 

X J2 I ^ } I [iW^ {si2S3)s] j) , (40) 

s 



where -j^ ^^j^ ^^J ^ is a 6-J symbol. 



For the states that we consider, the explicit expressions for the HQET states in terms of the quark model 
states are 

|l/2-,j = l) = y||l/2-,L = l,5-l/2) + -i=|l/2-,i=l,5 = 3/2), 
|l/2-,j = 0) = --^|l/2-,L=l,5=l/2) + y||l/2-,L = l,5 = 3/2), 
|3/2-,j-2) - y||3/2-,L = l,5=l/2) + i=|3/2-,L=l,5 = 3/2), 



|3/2-,j 



1) = --^|3/2-,L=l,5=l/2) + y||3/2-,L = l,5 = 3/2), 



|3/2+,j = 2) = (|3/2+,L = 2,5" = 1/2) + |3/2+,L = 2,5" 3/2)) , 

|3/2+,j = l) = i=(-|3/2+,L = 2,^=l/2) + |3/2+,L = 2,5 = 3/2)), 

|5/2+,j = 3) = ^|5/2+,L = 2,5 = l/2) + ^|5/2+,L = 2,5 = 3/2), 

|5/2+,j = 2) = -^|5/2+,L = 2,5=l/2) + ^|5/2+,L = 2,5 = 3/2). (41) 
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For all of the quark model states shown on the r.h.s of these equations, 5=1/2 corresponds to spin wave 
function of the type. The form factors that describe transitions to these states are shown in Appendix |0 
Other states not shown above are single component states in both representations, and these are 



|l/2+,J 


= 1) = 


|l/2+,L = 


0,5 = 


1/2), 


|3/2+,j 


= 1) = 


|3/2+,L = 


0,5 = 


3/2), 


|1/2^J 


= 1) = 


1/2+ L = 


0,5 = 


1/2), 


|3/2+,.? 


= 1) = 


|3/2+ L = 


0,5 = 


3/2), 


|l/2j,J 


= 1) = 


|l/22+,i = 


2,5 = 


3/2), 


|5/2-,j 


= 2) = 


|5/2-,L = 


1,5 = 


3/2). 



The subscripts '1' denotes the first radially-excited copy of the ground state multiplet. The '2' denotes an 
orbitally-excited state with = 1/2+ . This state forms a j = 1 multiplet with the second 3/2+ state listed in 
Eq. (HU. 

We now examine the form factors of Appendix [HI along with some of the form factors in Appendix ^ and 
compare these with the predictions of HQET shown in Section IlII Bl We begin with a discussion of the decays 
to pseudotensor final states. 



1/2- 



The HQET predictions for decays to this state are shown in Eq. Hll|l , while the quark model form factors are 
shown in Section IC II Noting that w — 1 w O (1/mg), the leading order predictions are that Fi = F2 = G2 = 0, 
and F3 = — G3 = Gi. The form factors of Section FC ll satisfv these relations, and allow us to identify 



(0) 

m 



in the harmonic oscillator models, or 



ax V alv 



exp 



3ml 



(0) 
Tl2 



^V2- 



5/2 



3 



in the Sturmian models. In these expressions, and in those that follow, 



and 



2 1 / 2 I 2 \ 



(ixx'=7i{(ix+(i'x)- 



(43) 

(44) 

(45) 
(46) 



B. (1/2+, 3/2+) 

The HQET predictions for decays to this pair of states are shown in Eqs. (|12() and (|13|l . The three multiplets 
with these quantum numbers are discussed separately. 



1. Ground State 

The quark model form factors for the ground state doublet are shown in Sections IB II and IB 81 Comparison 
of these form factors with the predictions of HQET leads to 

- -\ri'^' (47) 
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at the non-recoil point, and allows us to identify 



(1) /aAOA'V^^ / Sml p2 



in the harmonic oscillator models, or 



3/2 

(49) 



in the Sturmian models. In the heavy quark limit, both forms yield the expected normalization at the non- recoil 
point, namely ri[^\w = 1) = —1. It must be emphasized that the relationship between r/^^"* and 772^'' given above 
is one that arises only in the context of the quark model. In HQET, these two Isgur-Wise functions are a 
priori independent of each other. A more complete expression of the relationship between T^j^'' and r]^^ can be 
obtained by noting that G2 and 6*3 for the 1/2+ final state both vanish at leading order in the quark model. 
This leads to 

v['\w) = -{l + w)vi'\w), (50) 
valid at leading order in the heavy quark expansion. 

2. Radial Excitation 

The form factors for decays to the radially excited (1/2+, 3/2+) multiplet are shown in Sections IB 2l and lB 91 
Comparison of these form factors with the predictions of HQET again leads to 

Vi''=-lv['\ (51) 



and allows us to identify 



(1) ^ / QAay y^' ( iml 



in the harmonic oscillator models, or 

3/2 



^1 - ~A 7 



4 Pxf3x' 



1 Jr- 3 "^l 

2 ,nf Pl^, 



(53) 



in the Sturmian models. As with decays to the ground state multiplet, the full relationship between jy^^-* and 
772^^ can be deduced to be 

,yW(u;) = -(1 + u))r/f'(u)), (54) 
valid at leading order in the heavy quark expansion. 

3. Orbital Excitation 



The orbitally excited (1/2+, 3/2+) multiplet has a very different structure from either of the two multiplets 
discussed previously, and the form factors that are non- vanishing at leading order are different. For the ground 
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state and radially excited multiplets, Fi, F2, F3 and Gi are the non- vanishing form factors at leading order for 
the 1/2+ state, for instance, and this pattern is repeated with the radially excited multiplet [ignoring, for the 
moment, the fact that a\ — a\, « 0{l/mg)]. For the orbitally excited states, whose form factors are shown 
in Sections IB 31 and IC 31 the pattern is different, with G2 and G3 being the non- vanishing form factors for the 
1/2+ state. 

Comparing these quark model form factors with the leading order predictions of HQET allows us to deduce 
that 



= 0, 



7/2 



exp 



3mi 



2™L 



(55) 



in the harmonic oscillator basis, or 



0, 



9 mi 



7/2 



3 rn; 
2 m'f- 



(56) 



in the Sturmian basis. 



C. (3/2-, 5/2-) 



The HQET predictions for this multiplet are shown in Eqs. (|14(l and (|15(l . while the quark model predictions 
for these states are shown in Sections IB 121 and IC 61 Comparison of these two sets of equations yields 




in the Sturmian models. 



D. 



(5/2+, 7/2+) 



The HQET predictions for this multiplet are shown in Eqs. (|15|) and (|17|l . while the quark model predictions 
for the 5/2+ state are shown in Section fC 81 We have not calculated the form factors for the 7/2+ state in our 
models. Comparison of the HQET predictions with the results of the quark model calculation yields 



(3), 



3 ml f axax' 



y/2 al \ a 



AA' 



7/2 



exp 



3mi 



(59) 



in the harmonic oscillator models, or 



rj^{^\w) = -{l + w)ri^^\w) = -3V3 



(3), 



PI 



3 

2 mil 



7/2 



(60) 



14 



in the Sturmian models. 

We now turn to a discussion of the decays to daughter baryons having tensor hght diquark. We note, first 
that there exists a 1/2+, j — singlet state. At leading order, the form factors for decays to such a state vanish 
in HQET. In the quark model, such a state can be constructed, but the overlap of its wave function with that of 
the decaying parent baryon is zero to the approximation to which we work, and is strongly suppressed beyond 
that. Thus we do not have form factors for such a state. For the remaining decays of the tensor type, there is 
a single Isgur-Wise type form factor. 

E. (1/2-, 3/2-) 

The HQET predictions for this multiplet are shown in Eqs. (|19l) and (|20|l . while the quark model predictions 
are shown in Sections IC 21 and KT^ Comparison of these two sets yields 



2™L «AA' 



(61) 



in the harmonic oscillator models, or 



r(i) = a/3- 



/3a 



5/2 



(62) 



in the Sturmian models. 



F. (3/2+, 5/2+) 

The HQET predictions for this multiplet are shown in Eqs. (|21|) and (|22ll . while the quark model predictions 
are shown in Sections IC 41 and IC 71 Comparison of these two sets yields 



«A 



(63) 



in the harmonic oscillator models, or 



^(2) =3^2^. 



V ) 



5/2 



PI 



1 4_ 3 rn j p2 



(64) 



in the Sturmian models. 



VI. NUMERICAL RESULTS 



A. Model Parameters, Mass Spectra and Wave Functions 



In Section llV Bl we introduced the two Hamiltonians we diagonalize to obtain the baryon spectrum. They 
differ only in the form chosen for the kinetic portion, one of which is non-relativistic (NR), while the other 
is semi-relativistic (SR). In addition, we use two different expansion bases to obtain the wave functions: the 
harmonic oscillator (HO) basis, and the Sturmian (ST) basis. In the following, the four spectra we obtain will 
be denoted HONR, HOSR, STNR and STSR, in what should be an obvious notation. 



15 



TABLE I: Hamiltonian parameters obtained from the four fits. In the first column, HO refers to the harmonic oscillator 
basis, while ST refers to the Sturmian basis. In the same column, NR and SR indicate non-relativistic and semi-relativistic 
Hamiltonians, respectively. The form of these Hamiltonians is described in Section fiV Bl 



model 


rricr 


nis 


nic 


nit 


b 


acoui 


Qhyp 


Cqqq 






(GeV) (GeV) (GeV) (GeV) (GeV^) 






(GeV) 




HONR 


0.39 


0.63 


1.90 


5.30 


0.16 


0.21 


1.18 


-1.50 


0.73 


HOSR 


0.39 


0.55 


1.81 


5.26 


0.13 


0.15 


0.86 


-1.11 


0.70 


STNR 


0.41 


0.63 


1.90 


5.30 


0.13 


0.23 


0.34 


-1.40 




STSR 


0.42 


0.60 


1.83 


5.31 


0.14 


0.08 


0.25 


-1.40 





There are eight free parameters (nine in HO models) to be obtained for each spectrum: four quark masses 
{iriu ~ rud, rris, rric and nib), and four parameters of the potential (ahyp, "Coui, b and Cqqq). These eight 
parameters are determined from a 'variational diagonalization' of the Hamiltonian. The variational parameters 
are the wave function size parameters ap and a\ of Eq. (|28() . or Pp and /?a of Eq. H29|l . This variational 
diagonahzation is accompanied by a fit to the known spectrum. In this fit, the eight parameters mentioned 
before are varied. In addition, it is important to include any experimentally known information from semileptonic 
decay rates into the fit of these parameters. At present, this information is limited to the decay rate for 
A+ Ke^v. The rationale here is that the dynamics leading to the spectrum of states also play a crucial 
role in the semileptonic processes we are studying. By incorporating known semileptonic decay rates in our 
fit, we expect that predictions for as yet unmeasured rates will be more robust. The values we obtain for the 
Hamiltonian parameters are shown in Table 

In I, we presented the values obtained for the fit parameters. In the present manuscript, we have modified our 
variational diagonalization procedure somewhat, so it is appropriate for us to show the resulting values of the 
fit parameters again. The modification of this procedure is best explained with a concrete example, that of the 
1/2+ baryons. The wave functions for such states are expanded in terms of the seven basis states of Eq. (I25II . 
When the Hamiltonian is diagonalized in this basis, we obtain wave functions for seven states with = 1/2+. 
The variational part of the computation can be carried out by minimizing the energy of any of these seven 
states. In I, we used the energy of the lowest-lying state for this variation, and this led to the choices for a\ and 
ap (or /3a and /3p) reported there. In the present work, we no longer use the lowest-lying state for the variational 
calculation, but one of the excited states. This means that the values of the wave-function size parameters, as 
well as of the parameters in the Hamiltonian, in addition to the compositions of the states, are modified from 
their values in I, sometimes significantly so. In our current work we have chosen the third lowest-lying state for 
the variational procedure. 

One consequence of including the decay rate for A+ Ae+j/ in the fit is that the light (u and d) quark masses 
we obtain are consistently larger than conventional values. Other quark masses do not differ significantly from 
those in our previous work. We also obtain somewhat different values for acoui and cthyp, and the strength of 
the hyperfine interaction in the two harmonic oscillator models is larger than we reported in I. As a result, the 
hyperfine splittings in the baryon spectrum are now quite well reproduced in the HONR and HOSR models. 
We also note that the values of 5, the slope of the linear potential, remain similar to our previous results, which 
tend to be smaller than in most published studies of hadron spectra. However, recent work by Barnes, Godfrey 
and Swanson [Tlj reports a value of 0.14 for this parameter, obtained by fitting a similar Hamiltonian to the 
spectrum of heavy mesons. 

In the two harmonic oscillator models we fit an additional parameter k, which appears in the A+ — > Ae+z/ 
form factors. We have discussed the origin of k in I. This parameter was introduced in an ad hoc manner in 
the model of Isgur, Scora, Grinstein and Wise (ISGW) to take into account "relativistic effects". In the 

harmonic oscillator models, all form factors are proportional to the exponential factor 
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which ISGW modify to 

We include this parameter in our calculation of the HO form factors and rates in part because the work we 
present is done in the same spirit as the the work of ISGW, and such a parameter was found to be necessary 
in Ref. However, instead of choosing a particular value, as was done in Ref. 12], we treat k as a free 

parameter constrained to lie between 0.7 and 1.0. The values we obtain for k are shown in Table for both 
HONR and HOSR models. We note, however, that we do not include this parameter in the form factors for the 
decays of heavy Q,q baryons. This parameter is meant to mimic relativistic effects in the spectator quarks in 
the decaying baryon, and such effects are expected to be smaller for s quarks than they are for u and d quarks. 

The effect of the parameter k < 1 is to soften the form factors. It has been established that nonrelativistic or 
semi-relativistic quark models using an oscillator basis tend to underestimate the charge radius of light-quark 
systems such as the proton, and that some part of this underestimation can be attributed to relativistic effects 
in the evaluation of the electromagnetic current [T^ . A procedure similar to the inclusion of the parameter k 
by ISGW was used by Foster and Hughes to modify electromagnetic form factors of light-quark systems 
calculated in a nonrelativistic quark model. 



TABLE II: Wave function size parameters, Op and a\, for states of selected with spin-flavor symmetric light diquark, 
in different models. All values are in GeV. 





model 








n 






(QA, Op) ("A, 


Op) 


(qa. 


Op) 


(oA, ctp) 


1/2+ 


HONR 


(0.54, 0.42) (0.49, 


0.42) 


(0.39, 


0.42) 


0.42 


1/2+ 


HOSR 


(0.51, 0.48) (0.49, 


0.47) 


(0.42, 


0.44) 


0.42 


1/2+ 


STNR 


(0.72, 0.42) (0.67, 


0.47) 


(0.49, 


0.65) 




1/2+ 


STSR 


(0.72, 0.38) (0.66, 


0.44) 


(0.63, 


0.78) 




3/2+ 


HONR 


(0.50, 


0.42) 


(0.37, 


0.40) 


0.40 


3/2+ 


HOSR 


(0.48, 


0.47) 


(0.36, 


0.41) 


0.40 


3/2+ 


STNR 


(0.74, 


0.40) 


(0.43, 


0.71) 




3/2+ 


STSR 


(0.67, 


0.43) 


(0.63, 


0.78) 




1/2- 


HONR 


(0.51, 


0.43) 


(0.38, 


0.41) 


0.42 


1/2- 


HOSR 


(0.48, 


0.47) 


(0.40, 


0.42) 


0.42 


1/2- 


STNR 


(0.67, 


0.47) 


(0.61, 


0.53) 




1/2- 


STSR 


(0.65, 


0.45) 


(0.60, 


0.75) 




5/2" 


HONR 


(0.51, 


0.43) 


(0.37, 


0.42) 


0.42 


5/2" 


HOSR 


(0.48, 


0.47) 


(0.39, 


0.42) 


0.42 


5/2- 


STNR 


(0.67, 


0.47) 


(0.61, 


0.53) 




5/2- 


STSR 


(0.65, 


0.45) 


(0.60, 


0.75) 





In carrying out our fits, we generally allow the values of ap to be different from ax, as in I. The exceptions 
occur in cases when the three quarks are identical, as they are in the nucleon or the SI. In such cases, the 
variational diagonalization automatically selects ap — a\ in the HO bases. Table Ull shows some of the values 
we obtain for the size parameters. The omitted parameters for the states that are significant for this work are 
related to those presented. For instance, for the 1/2^ states, the size parameters are the same as for the 1/2+ 
states. Furthermore, since we do not include a spin-orbit interaction in our Hamiltonian, the size parameters for 
the 1/2- and 3/2- states are identical. We do not show the size parameters for the S states with = 5/2+ 
or 7/2+ mainly because we find that semileptonic decay rates to these states are very small. We also omit the 
size parameters for the analogous Qq states with Q = c, s. 
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B. Mass Spectra 

Portions of the mass spectra we obtain using our four models are shown in Tables IIIII and IIVI In these 
tables, the first two columns identify the state and its experimental mass, while the next four columns show the 
model masses that result from a fit of the Hamiltonian parameters to those states whose experimental masses 
are known. We note that for the Q and S states, the predicted masses are in satisfactory agreement with the 
available experimental values, with little variation among the results from the different models for these states. 



TABLE III: Baryon masses in GeV in the quark models we use. Hamiltonian parameters for each model are obtained 
from fits to the experimental masses where known; other masses shown are predictions of the models. The first two 
columns identify the state and its experimental mass, while the next four columns show the masses that result from the 
models. 



State 


Experimental Mass 


HONR HOSR STNR STSR 


H(l/2+) 


1.32 


1.32 


1.35 


1.40 


1.39 


H(l/2+(rad)) 




2.03 


1.79 


2.06 


1.83 


H(l/2+(or6)) 




2.08 


2.00 


2.10 


1.95 


H(3/2+) 


1.53 


1.52 


1.54 


1.45 


1.46 


H(3/2+(rad)) 




2.16 


2.08 


2.10 


2.12 


H(3/2+(or6)) 




2.14 


2.18 


1.98 


1.96 


H(3/2-) 


1.82 


1.83 


1.78 


1.79 


1.80 


H(5/2-) 




1.84 


1.78 


1.78 


1.81 


H(5/2+) 




2.08 


2.14 


2.08 


2.00 


n(3/2+) 


1.67 


1.66 


1.66 


1.60 


1.67 


Q{3/2+ (rad)) 




2.20 


2.07 


2.34 


2.13 


Q{3/2+ (orb)) 




2.23 


2.11 


2.24 


2.14 


f2(3/2-) 




1.95 


1.84 


1.88 


1.88 


fl{5/2~) 




1.95 


1.89 


1.89 


1.89 


nc(l/2+) 


2.70 


2.69 


2.72 


2.73 


2.71 


nc{l/2+ (rad)) 




3.18 


3.09 


3.24 


3.24 


flc{l/2+{orb)) 




3.25 


3.17 


3.24 


3.26 


nc(3/2+) 




2.77 


2.78 


2.75 


2.73 


n^{3/2+{rad)) 




3.22 


3.15 


3.30 


3.24 


nc{3/2+ (orb)) 




3.24 


3.18 


3.23 


3.26 


n,{i/2-) 


3.00 


3.00 


2.97 


3.00 


3.02 


^^(5/2-) 




3.02 


2.99 


3.01 


3.02 


!^i,(l/2+) 




6.08 


6.13 


6.08 


6.14 



In Table lrvl we also present some of the masses of the nucleons and Ag states, mainly to show the improvement 
that has resulted from the modified variational procedure. We have obtained a better spectrum for almost all of 
the nucleons and Ag baryons, with significant improvement in the iV(1440) and the A resonance model masses. 



C. Wave Functions 

Significant mixing of wave function components occurs in many of the Qq and S states, for all flavors, partic- 
ularly in the Sturmian models. The mixing coefficients that result, along with recalculated mixing coefficients 
for N and Aq states, are tabulated in Tables El and IVII for all four models. In Table we show the wave 
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TABLE IV: Baryon masses in GeV fitted in the four quark modeis we use. Tlie first two coiumns identify tlie state and 
its experimental mass, wfiile tfie next four columns show the masses that result from the models. 



State 


Experimental Mass 


HONR HOSR STNR STSR 


N{l/2+) 


0.94 


1.00 


1.08 


1.07 


1.12 


N{l/2+) 


1.44 


1.68 


1.56 


1.76 


1.58 


N{l/2-) 


1.54 


1.47 


1.47 


1.51 


1.47 


7V(3/2+) 


1.72 


1.72 


1.76 


1.77 


1.73 


A(3/2+) 


1.23 


1.24 


1.32 


1.20 


1.20 


A(l/2+) 


1.12 


1.11 


1.11 


1.09 


1.05 


A(l/2+) 


1.60 


1.74 


1.63 


1.61 


1.59 


A(l/2-) 


1.41 


1.49 


1.50 


1.46 


1.52 


A(3/2+) 


1.89 


1.85 


1.74 


1.73 


1.81 


Ac(l/2+) 


2.28 


2.27 


2.26 


2.27 


2.21 


A,(l/2-) 


2.59 


2.63 


2.60 


2.60 


2.66 


A,(l/2+) 


5.62 


5.62 


5.62 


5.62 


5.62 



function coefficients for the 1/2+, 3/2+ and 3/2~ states, in each flavor sector, for each model for the VIq and S 
baryons. The exceptions are for O, where we do not use the Sturmian basis. We do not present the mixing for 
1/2^ states because they have exactly the same mixing coefficients as, and are degenerate with, the 3/2~ states. 
For other states we treat, such as l/2|j~, 3/2|j~ and 5/2~, the wave functions that result are single component 
wave functions. The mixing shown in these tables complicates the extraction of the form factors. However, in 
all numerical results that we show for the form factors and the decay rates, this mixing is properly taken into 
account. 



TABLE V: Mixing coefficients (77^) of the lowest-lying 1/2+, 3/2+ and 3/2 states of and Q,q, in different fiavor sectors. 
The rii axe defined in Appendix 1X1 



Baryon states 




HONR 






HOSR 




STNR 






STSR 








r)i 


ri2 


r)3 


r?i 


772 




r)i 


J72 


7?3 






»73 


2(1/2+) 


0.970 


0.100 


0.198 


0.962 


0.062 


0.230 


0.969 


-0.226 


0.093 


0.964 


-0.256 


0.058 


S(3/2+) 


0.996 


0.077 


0.033 


0.999 


-0.009 


0.038 


0.947 


-0.313 


0.066 


0.935 


-0.334 


-0.118 


H(3/2-) 


0.484 




0.875 


0.641 




0.767 


-0.296 




0.955 


0.115 




0.993 


n,(l/2+) 


0.976 


0.093 


0.189 


0.980 


-0.035 


0.189 


0.980 


0.200 


0.025 


0.933 


0.361 


\r)3\ < 0.001 


nc(3/2+) 


0.995 


0.061 


0.072 


0.993 


-0.091 


0.068 


0.964 


0.243 


-0.107 


0.948 


0.2317 


-0.010 


nc(3/2-) 


-0.234 




0.997 


-0.293 




0.956 






1.00 






1.00 


nb(l/2+) 


0.985 


0.086 


0.147 


0.980 


-0.09 


0.173 


0.957 


0.291 


0.006 


0.937 


0.350 


0.005 






HONR 






HOSR 


















J?! 






J?! 


m 


n-i 














0(3/2+) 


0.996 


0.063 


0.063 


0.998 


0.042 


0.042 
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TABLE VI: Mixing coefficients (77^) of tlie lowest-lying 1/2^ states of N and Ag in different flavor sectors. The rji are 
defined in Appendix IXI 



Baryon states 


HONR 

vi m 


HOSR 
rii J72 Vs 


STNR 

7?1 V2 V'i 


STSR 
Vi V2 V:i 


N{l/2+) 
A(l/2+) 
A,(l/2+) 
A6(l/2+) 


0.959 0.095 0.246 
0.976 0.186 0.026 
0.976 0.186 0.103 
0.977 0.185 0.106 


0.949 0.067 0.272 
0.933 0.289 -0.089 
0.939 0.337 0.006 
0.929 0.368 0.024 


0.963 0.219 0.154 
0.962 0.224 0.158 
0.951 0.227 0.206 


0.869 0.478 0.129 
0.876 0.447 0.179 
0.861 0.438 0.259 



D. Form Factors and Decay Rates: CIq 

1. VLb ni*^ 

The form factors at the non-recoil point for the flf, ^ flc decays calculated using all four models are shown 
in Table IVTll We show only form factors for the decays to final states with = 1/2+ and 3/2+. These states 
constitute the elastic channels. We note that the form factor values at the non-recoil points are very close to 
each other in all our models. 

It is instructive to compare our form factor values at the non-recoil point with the HQET predictions for 
decays to the ground state doublet. For the state with = 1/2+, the HQET prediction is Fi = Gi — Jyi^V^- 
With the known normalization condition on ri[^\ this means that we expect Fi — Gi = —1/3 at the non-recoil 
point. Our results for Gi are very close to this, but those for Fi are not. The deviations from the HQET 
predictions can easily be traced to the presence of l/m^ and l/mg terms in the quark model results for i^i, 
while no such terms exist in the model predictions for Gi. If such terms are ignored in Fi then the HQET 
prediction is indeed satisfied in the quark models. Similarly, HQET predicts that F2 and F3 should each have 
the value of 2/3 at the non-recoil point. The model results agree with this prediction for F2 but not for F-^, 
and the differences can again be traced to the presence of 1/m^ terms in the quark model results. The HQET 
predictions for G2 and G3 are less easily interpreted, but comparison of those predictions with the quark model 
calculations suggests that the Isgur-Wise function 772^^ is normalized to 1/2 at the non-recoil point. We have 
already raised this point in Section IVI where we compare the quark model predictions for the form factors with 
those of HQET. 

For the state with = 3/2+, HQET predicts that F4 = -G4 = 2/\/3. This is well satisfied by the model 
predictions for G4, but the model predictions for include l/rriq contributions. Assuming that 772^^ is indeed 
normalized to 1/2 at the non-recoil point, the HQET prediction is then that Gi = and F3 = — G3 = — I/a/S- 
The quark model results for G3 are close to the HQET prediction, but those for F3 deviate from this prediction 
because of the presence of l/rriq terms in the quark model results. 



TABLE VII: The form factors for Qb transitions to the ground-state Qc multiplet with — (1/2+, 3/2+), calculated 
at the non-recoil point, in the four models used here. 



spin 


model 


Fi F2 


F3 


F4 


Gi G2 


G3 


G4 


1/2+ 


HONR 


-0.48 0.58 


0.86 




-0.32 0.13 


-0.03 




1/2+ 


HOSR 


-0.47 0.56 


0.87 




-0.32 0.12 


-0.02 




1/2+ 


STNR 


-0.47 0.63 


0.83 




-0.33 0.11 


-0.04 




1/2+ 


STSR 


-0.47 0.62 


0.83 




-0.33 0.11 


-0.04 




3/2+ 


HONR 


0.80 0.0 


-0.80 


1.61 


0.0 -0.17 0.63 


-1.12 


3/2+ 


HOSR 


0.78 0.0 


-0.78 


1.58 


0.0 -0.16 


0.62 


-1.12 


3/2+ 


STNR 


0.83 0.0 


-0.83 


1.65 


0.0 -0.18 


0.64 


-1.15 


3/2+ 


STSR 


0.82 0.0 


-0.82 


1.64 


0.0 -0.19 


0.63 


-1.14 
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FIG. 1: Form factors for Vih Qc{l/2'^) obtained using harmonic oscillator wave functions (left panel, HOSR and 
HONR models) and Sturmian wave functions (right panel, STSR and STSR models). In each panel, the solid curves 
arise from the semi-relativistic version of the model, while the dashed curves arise from the non-relativistic version. 



Figure^shows the dependence of the form factors for the elastic transition flj, ilc{l/2^) calculated in 
the HONR and HOSR models on the left, and in the STSR and STNR models on the right. In each panel, the 
solid curves arise from the SR version of the model, while the dashed curves are from the NR version. Here we 
note that the form factors calculated using the Sturmian wave functions have slopes near the non-recoil point 
that are similar to those calculated using the harmonic oscillator wave functions. This is due to the fact that 
we have similar mixing patterns for the Qb and Qc ground state wave functions in all models as can be seen in 
Table El 

The differential decay rates dT/dq^ obtained in the four models for different final states in fif, Vt^c^ Iv^, with 
i = e~,/i~, are shown in the upper panels of Figure[21 For these rates, we use \Vcb\ = 0.041. In these figures, 
we only show the differential rates for the dominant decays to the two elastic channels, with = 3/2+ and 
— 1/2+, and for two orbital excitations, the states with = 3/2^ and 5/2^. We have also examined the 
differential decay rates to the l/2~, l/2j~, and 3/2^ orbitally excited states, as well as to the radially excited 
states 1/2]'" and 3/2^'" (notations defined in Section FlV All . We have found that the branching fraction for the 
radially excited states (not shown in Table rVlTT|l are small, whereas the branching fraction for the decays to 
the orbitally excited states are not insignificant, as shown in Table IVlTTl The lower panels of Figure |21 show the 
differential decay rates of f2f, decaying to the same Vic final states as in the upper panels, but with a r lepton 
in the final state. 

In Table rvTTTl we show the integrated decay rates obtained for the selected final states in the four quark models 
we use. The first part of this table shows the rate with a vanishing lepton mass, while the second part shows 
the rate when the final lepton is a r. The last two rows of the first part of the table present the total decay 
rate and the ratio of the elastic to the total semileptonic decay rate. The integrated rates for the elastic decay 
modes (l/2+,3/2+) obtained in all models are similar. However, the two Sturmian models predict somewhat 
smaller rates for decays into the inelastic flc channels. As a result the branching fraction for the elastic decay 
mode is smaller in the HO models than in the Sturmian models. If we consider the two HO models alone, the 
predicted elastic branching fraction is 49.5±1.5%. The corresponding prediction from the Sturmian models is 
67.5±0.5%. Thus, the two HO models are consistent with each other, and the two ST models are consistent 
with each other, but the HO and ST models are in disagreement. Both sets of models predict that the elastic 
decay processes dominate the fit semileptonic decay but do not saturate it; there is some significant branching 
fraction to the inelastic channels. 
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FIG. 2: The differential decay rates for selected Qb transitions, in the various models that we use. The panels 

on the left arise from the two versions of the harmonic oscillator model, while those on the right are from the Sturmian 
models. The upper panels are for Qt ^li'^t^ui, where £ is e~ or fi~ . The lower panels are for Qi, — > r2i*Vi/T. The curves 
are for final states with = 3/2'^ , 1/2'^ , 3/2~ and 5/2~. In each panel, the solid curves arise from the semi-relativistic 
version of the model, while the dashed curves arise from the non-relativistic version. 



We may also compare our decay rates with HQET predictions. As discussed in Section Fill Bl there are a 
number of pairs of degenerate states, such as states with = (l/2+,3/2+), = (1/2^,3/2^) and = 
(3/2j~, 5/2^). In the heavy quark limit, the ratio between the rates of the ground state heavy baryon decaying 
to the states in the first two degenerate pairs is expected to be 1 : 2, and for the third pair it is 2 : 3. In 
other words, for example, we expect the rate for $1^(1/2+) ric(3/2+) to be twice as large as the rate for 
r2h(l/2+) $7^(1/2'''). The expected pattern of the rates coming from the HQET prediction is reflected in all 
of our quark model calculations, as can be seen in Table IVIIII Departures from these predictions are due to 
l/mg and l/niq corrections, and the fact that, for instance, the 3/2~ state shown in the table is not exactly 
the state in the (l/2~,3/2^) multiplet, but contains some admixture of the 3/2~ state from the (3/2j~,5/2~) 
multiplct. 

In Table IVIlTl we have shown only the rates to the fli*^ states which have a significant branching fraction. 
However, we have calculated the rates of fib decaying to the majority of states listed in Eq.|^ and have found 
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TABLE VIII: Integrated decay rates for Qb f^c*' in units of lO^^s^^, for a selection of flc states in the four models 
we consider. The upper portion of the table presents the decay rates obtained with massless leptons (electron or muon), 
while the lower portion corresponds to the rates with a massive r lepton in the final state. 





fib ^ O.i'h-pe 


Spin 


r(HONR) r(HOSR) 


r(STNR) r(STSR) 


3/2+ 


1.68 


1.93 


2.01 


2.00 


1/2+ 


0.71 


0.94 


0.87 


0.96 


3/2" 


0.69 


0.71 


0.43 


0.43 


3/2r 


0.44 


0.47 


0.23 


0.22 


5/2" 


0.70 


0.72 


0.46 


0.45 


1/2- 


0.32 


0.33 


0.20 


0.20 


1/2- 


0.44 


0.48 


0.07 


0.07 


total 


4.98 


5.58 


4.27 


4.33 


rn^(l/2+,3/2+)/rtotal 


0.48 


0.51 


0.67 


0.68 




fit ni*V T^T 


Spin 


r(HONR) r(HOSR) 


r(STNR) r(STSR) 


3/2+ 


0.81 


0.89 


0.89 


0.88 


1/2+ 


0.42 


0.49 


0.45 


0.49 


3/2- 


0.15 


0.14 


0.10 


0.11 


3/2r 


0.09 


0.10 


0.07 


0.07 


5/2- 


0.14 


0.14 


0.11 


0.10 


1/2- 


0.07 


0.06 


0.05 


0.05 


1/2- 


0.08 


0.09 


0.02 


0.02 


total 


1.76 


1.91 


1.69 


1.72 



that the rates for states not shown in Table IVTllI arc small (of the order of 1% of the total rate that we have 
obtained). 

The differential decay rates dT/dq^ for flc decaying to the dominant final states of S*^*) are shown as functions 
of in Figure|21 Here, we use \Vcs\ = 0.224. The left panel shows the differential decay rates obtained in the 
HONR and HOSR models, and the right panel shows the results from the STNR and STSR models. In both 
harmonic oscillator and Sturmian models we see significant differences between the non-relativistic and semi- 
relativistic predictions, with the most significant difference occuring between the HOSR and HONR predictions 
for the 1/2+ final state. Although the different size parameters obtained in various models play a role in these 
differences, we should also note here that this particular decay rate is very sensitive to the value of rria {a = s 
in this decay). Our fitted values for the strange quark mass in the different models, shown in Tabled show 
significant variation. We have seen similar differences between NR and SR model differential decay rates for 
the Ac ^ n decay mode presented in I. 



The integrated decay rates for flc decaying to various states of S are shown in Table IIXI in the four different 
models we use. The last two rows give the total decay rate and the branching fraction to the 'elastic decay' 
channels. The two elastic decay modes dominate the decay, with a small fraction of ilc decays going to the 
orbitally excited S^*^ states. The rates obtained in the Sturmian models are smaller than those obtained in the 
harmonic oscillator models. Nevertheless, the predicted elastic fraction does not depend strongly on the models 
we use. 
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FIG. 3: The differential decay rates for selected f2c — > H'*-* transitions, in the various models that we use. The panel on 
the left arises from the harmonic oscillator models, while the one on the right arises from the Sturmian models. They 
key to the curves is as in Fig. |5| 

TABLE IX: Integrated decay rates for S7c H'*' in units of 10^" s~^, for selected H states in the four models we consider. 



Spin 


r(HONR) r(HOSR) 


r(STNR) r(STSR) 


3/2+ 


0.51 


0.59 


0.36 


0.39 


1/2+ 


0.39 


0.65 


0.34 


0.39 


3/2" 


0.05 


0.05 


0.03 


0.02 


5/2" 


0.03 


0.03 


0.04 


0.03 


1/2- 


0.04 


0.04 


0.05 


0.05 


1/2- 


0.08 


0.11 


0.02 


0.02 


total 


1.10 


1.47 


0.84 


0.90 


rH(l/2+,3/2+)/rtotal 


0.82 


0.84 


0.83 


0.87 



In Table Ixl the integrated decay rates of fib to various states of S are presented, along with the rates for the 
rjfc — > 'E!^*'^t^ Vr in the second part of the table. We note that due to the large phase space available, a significant 
fraction of f2f, is predicted to decay semileptonically to various excited S states. Apart from a few decay modes, 
such as E!(1/2J"), the integrated decay rates of Jib to the various states of S vary minimally within the models 
we use. 



3. Qc 



Table IXll shows the integrated rates for decaying to in our two harmonic oscillator models. The decay 
to the ground state {J^ = 3/2+) almost saturates the decay, providing 91% of the total decay rate. The lack 
of phase space suppresses the decays to 'inelastic' channels. 
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TABLE X: Integrated decay rates for Slj, H^*' in units of 10^'^|K,6ps^^, for a selection of H states in the four models we 
consider. The last row of the top portion of the table shows the 'elastic fraction' (evaluated using the lowest lying 1/2+ 
and 3/2"^ states) obtained in our model, where the decays shown in the table are assumed to saturate the semileptonic 
decays. 













Spin 


r(HONR) r(HOSR) 


r(STNR) r(STSR) 


3/2+ 


2.90 


3.74 


2.16 


3.60 


1/2+ 


0.82 


1.46 


1.09 


1.78 


3/2" 


2.17 


2.73 


3.80 


3.96 


3/2- 


2.25 


4.15 


2.04 


2.03 


5/2" 


2.15 


3.73 


3.27 


3.20 


1/2- 


0.64 


0.87 


1.33 


1.39 


i/2r 


1.15 


1.80 


0.35 


0.37 


total 


12.08 


18.48 


14.04 


16.33 


rH(l/2+,3/2+)/rtotal 


0.31 


0.28 


0.23 


0.33 












3/2+ 


2.42 


3.10 


1.76 


2.83 


1/2+ 


0.89 


1.53 


1.06 


1.61 


3/2" 


1.86 


2.32 


2.88 


2.99 


3/2r 


1.86 


3.36 


1.39 


1.35 


5/2" 


1.62 


2.75 


2.33 


2.25 


1/2- 


0.59 


0.77 


1.04 


1.07 


1/2- 


1.09 


1.67 


0.30 


0.30 


Total 


10.33 


15.50 


10.76 


12.4 



TABLE XL Integrated decay rates for Qc Si'*-' in units of 10^^ s"^, for a selection of D, states in the two HO models 
we consider. The last row shows the 'elastic fraction' obtained in our model, assuming the decays shown in the table 
saturate the semileptonic decays. 



Spin 


r(HONR) r(HOSR) 


3/2+ 


4.74 


5.39 


3/2" 


0.31 


0.33 


1/2- 


0.18 


0.19 


total 


5.23 


5.91 


rn{3/2+)/rtotal 


0.91 


0.91 



E. Form Factors and Decay Rates: Aq, revisited 

One consequence of the modified variational procedure mentioned in Section IVI Al coupled with fitting to 
the decay rate for A+ Ae+i^, is that most of the decay rates for Aq have been modified. In this section, we 
briefly discuss the new results for Aq decays, and compare them with those presented in I. We also discuss the 
effects of K on these decay rates. 



Table IXU shows the integrated decay rates for Ac A^*^ obtained in the various models we use. The second 
and third columns of this table show the rates we obtain in the two harmonic oscillator models. Note that the 
rates obtained in the HO models are somewhat improved from those presented in I, but are still larger than the 
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TABLE XII: Integrated decay rates for Ac A'*' in units of lO^^s^^, in the models we consider. The row labeled 'total' 
is obtained by adding all calculated exclusive decay rates. The fourth and fifth columns show the K-affected integrated 
rates in the HO models. Integrated rates from our previous work are also shown in parentheses for comparison. The 
last row shows the 'elastic fraction' estimated in our models, assuming that the exclusive channels calculated saturate 
the semileptonic decays of the Ac. In each row, the numbers in parentheses are from I 





r(HONR) 


r(HOSR) 


r(HONRK) 


r(HOSRK) 


r(STNR) 


r(STSR) 


Expt. fl6] 


1/2+ 


1.80 (2.10) 


2.16 (2.36) 


1.07 


1.31 


1.32 (0.79) 


1.44 (1.11) 


1.05±0.35 


total 


2.00 (2.36) 


2.37 (2.73) 


1.22 


1.47 


1.53 (0.97) 


1.60 (1.31) 




TA/Ttotal 


0.90 (0.89) 


0.91 (0.86) 


0.88 


0.89 


0.86 (0.81) 


0.90 (0.85) 


1.0 (assumed) 



experimentally measured rate by about a factor of two, even though the rate is included in the fit. The fourth 
and fifth columns show the rates we obtain when k is included in the harmonic oscillator form factors. The 
second row in the table shows the rates we obtain with the modified variational procedure, while the numbers 
in parentheses in that row are the corresponding rates presented in I. In the Sturmian models, the elastic rate is 
larger than in I, resulting in elastic fractions that are similar to those obtained in the HO models in the present 
work. As a result, the overall model dependence of the elastic fraction has decreased significantly. One of the 
reasons for this improvement is that with the new minimization scheme both harmonic oscillator and Sturmian 
model wave functions for different baryon states have similar mixing coefficients, as can be seen in Tables fVl 
and El 

The third row presents the total rate we estimate for A+ A^*'>e~^i' (the numbers in parentheses are from I). 
Examination of these rates indicates that the model dependence in the A+ Ae+z/ rate has decreased compared 
with I. We have obtained an average value of 1.28]']q'2i x lO^^s^^ for the A+ Ke^v rate in the Sturmian models 
and the two K-modified HO models, consistent with the value reported by the CLEO collaboration H^. In the 
last row of Table IXIII the 'elastic fractions' obtained in the various models are shown. Our new calculations 
predict smaller inelastic branching fractions (9 to 14%) than reported in I (11 to 19%). The elastic fraction, 
averaged over the Sturmian models and the K-modified harmonic-oscillator models, is found to be 0.88 ± 0.02. 



TABLE XIII: Rates for At Ac*' decays in units of lO^^s"^. The fourth and fifth columns show the K-affected 
integrated rates in the HO models. Integrated rates from our previous work are also shown in parentheses for comparison. 
The row labeled 'total' is obtained by adding all calculated exclusive decay rates, while the row with the branching 
fractions assumes that the exclusive channels calculated saturate the semileptonic decays of the At. The elastic fraction 
reported by the DELPHI collaboration 17] (third row of numbers, eighth column) is actually r(At, — > Kctvt) /\y{^b — > 
Kclvi) + r(A(, — » AcTTTr^S-f)]. The errors on both numbers from DELPHI are statistical and systematic, respectively. 





r(HONR) 


r(HOSR) 


r(HONRK) 


r(HOSRK) 


r(STNR) 


r(STSR) 


Tdelphi [17] 


1/2+ 


3.50 (4.60) 


4.14 (5.39) 


1.81 


2.29 


2.39 (1.47) 


2.52 (2.00) 


/, r|7+0-90+1.30 
-0.65-0.98 


Total (Ai*'^~i?f) 


4.83 (5.95) 


5.21 (6.82) 


2.49 


2.89 


3.31 (2.36) 


3.04 (2.88) 




TAc /Ttotal 


0.72 (0.76) 


0.79 (0.79) 


0.73 


0.79 


0.72 (0.62) 


0.83 (0.69) 


r, ,,7+0.10+0.07 
-0.08-0.06 



The integrated rates for At Ac are shown in Table IXlTll Here we again present both the results from this 
work and from I for comparison. We note that the HO and ST models now provide predictions for both the 
elastic Af, — > Ac decay and the total decay rates that are more consistent with each other. It is intriguing that, 
for the decays of both the A;, and Ac, the results obtained in the K-modified HO models are consistent with the 
results obtained in the ST models. 

The elastic fraction for the A;, ^ Ac decay obtained in the various models is shown in the last row of 
Table IXIIII In the present work this fraction is slightly less model-dependent than in I. The average value 
for the Afc — > Ace~z^ rate obtained using the two Sturmian models and the two K-modified HO models is 
2.251 

a44 ^ 10^°s ^. These four models predict an average value of 0.77 ± 0.05 for the elastic fraction. 
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TABLE XIV: Decay rates of At N'^*^+rue in units of 10^^ |K,i,ps"^ The fourth, fifth, eighth and ninth columns 
show the K-affected integrated rates in the HO models. The row labeled 'total' is obtained by summing all the exclusive 
decay rates shown in the table. Integrated rates from our previous work are also shown in parentheses for comparison. 
The first four columns of numbers are for decays with a muon or electron in the final state, while the last four columns 
are for decays with a r in the final state. Also shown are the rates for Ac — > N^*^'^l^i'e in units of 10^°s~'^. 





Ab N^*'>+rve 


At N'-*'>+T-i?r 




r(HONR) 


r(HOSR) 


r(HONRK) 


r(HOSRK) 


r(HONR) 


r(HOSR) 


r(HONRK) 


r(HOSRK) 


1/2+ 


2.38 (4.55) 


3.53 (7.55) 


0.90 


1.19 


2.12 (4.01) 


3.11 (6.55) 


0.81 


1.06 


1/2+ 


0.44 


0.85 


0.49 


0.65 


0.35 


0.64 


0.12 


0.22 


1/2- 


1.18 


1.14 


0.45 


0.42 


0.95 


0.85 


0.38 


0.35 


3/2" 


0.87 


0.73 


0.20 


0.15 


0.64 


0.46 


0.17 


0.06 


3/2+ 


0.52 


0.28 


0.11 


0.06 


0.36 


0.17 


0.09 


0.05 


5/2+ 


1.29 


0.44 


0.18 


0.06 


0.69 


0.18 


0.13 


0.04 


Total 


6.68 (12.25) 


6.97 (21.31) 


2.33 


2.53 


5.11 (9.00) 


5.41 (15.53) 


1.70 


1.78 




Ac N^*^°£+ue 










1/2+ 


0.68 


0.84 


0.31 


0.38 










1/2- 


0.03 


0.03 


0.02 


0.02 











In Table IXTVl we present the integrated rates for Af, N^*''^£~'De obtained in the two HO models. Integrated 
rates for Af, — > N'^*^^t~ i>r are also shown in this table. From the integrated rates for all allowed decay modes to 
ground state nucleons and N* states shown in Table IXTVl we predict that a significant fraction of these decays 
are to excited nucleons, due in part to the large phase space available. We also present ^ p rates and total 
decay rates from I in parentheses for comparison. We note that, as in the two previous decay processes, the 
modified variational procedure leads to results in the two HO models that are more consistent with each other. 

VII. CONCLUSION AND OUTLOOK 

A constituent quark model calculation of semileptonic decays of flQ, which has several intriguing features, has 
been described in the preceding sections. Analytic results for the form factors for the decays to ground states 
and excited states with a selection of quantum numbers are evaluated, and compared to HQET predictions. 
For J7b rtc transitions, the HQET relationships among the form factors detailed in Appendix [21 in the 
(1/2+, 3/2+), (1/2-, 3/2"), (3/2", 5/2"), and (3/2+, 5/2+) doublets are respected in our predicted form 
factors, at leading order. 

These form factors depend on the size parameters of the initial and final baryon wave functions, and so a fit is 
made to the spectrum of the states treated here. Two model Hamiltonians are used, with either a nonrelativistic 
or semi-relativistic kinetic energy term, and with Coulomb and spin-spin contact interactions, as we did before 
in our work on Ag decay. For the present work we have modified our variational procedure from that used in I, 
and we have incorporated the A+ ~* Ae+z^ rate in the fit to the spectrum. As a result our size parameters are 
somewhat different from those shown in I, and our four baryon spectra are improved from those shown in I. 

As in I, the wave functions are expanded in either a harmonic oscillator or Sturmian basis, up to second-order 
polynomials, and our numerical results for form factors and rates are calculated using the resulting mixed wave 
functions. Four sets of predictions are made for form factors and rates, with wave functions, size parameters and 
mixing coefficients arising from fits using both the non-relativistic and semi-relativistic Hamiltonians, and using 
the two different bases. The variation among these predictions can be used to assess the model dependence in 
the results we obtain. 

At present there exist few quantitative measurements of the semileptonic decay rates of fig. The CLEO 
collaboration has measured B{n° n^e+i^) ■ a{e+e- n^X) = 42.2 ± 14.1 ± 11.9 fb. In addition the 
ARGUS and BELLE |8|| collaborations have reported the fl° (}~ decay as 'seen'. However, the branching 
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fraction for this decay has not yet been measured, which means that a comparison of our model predictions with 
experiment is not possible at this time. Nevertheless, it is instructive to examine our predictions for the various 
integrated rates for fit, and Oc decays in light of HQET predictions for these decays. There are a number of pairs 
of degenerate states, such as states with = (1/2+, 3/2+), = (1/2", 3/2"), = (3/2;,", 5/2"), for which 
HQET predicts the relationships between the form factors as well as the ratio of rates of the degenerate pairs. 
According to this prediction, the rate for the fit — > i7c(3/2+) is expected to be twice as large as the rate for 
fib — > fic(l/2+), and the rate for the = 3/2~ final state is also expected to be twice as large as the rate for 
the = l/2~ final state decay. In the same way the ratio of the rates for the final states of the (3/2j~, 5/2~) 
doublet is expected to be 2:3. The states we obtained in our spectral fits are not exactly the HQET states, but 
appear to be close approximations to such states, because these relationships among the rates of the degenerate 
pairs are well respected in all our quark model calculations for the semileptonic decays of fib- 

Our predictions for the semileptonic elastic branching fraction of fig vary minimally within the models we 
use. We obtain an average value of (84± 2%) for the elastic fraction of flc — > S^*'. For the case of Qc 
91% of all decays are elastic. The two HO models give results that are consistent with each other (48%, 51%) 

for the elastic fraction of £1^ — > . The two Sturmian models also give results that are consistent with each 
other (67% and 68%), but which are somewhat different from the fraction predicted by the HO models. 

We have recalculated the Aq rates in all four models. Following ISGW, we have also modified the results of 
the HO models by including a parameter k in the exponential factor that multiplies all HO form factors. As a 
result the model dependence of the Aq rates and the elastic rates have been decreased significantly. We have 
obtained an average value of 1.28lg'2i ^ lO^^s^^ for the A+ Ke^v decay rate in the Sturmian and K-modified 
HO models, and the elastic fraction for the decays of the Ac is 0.88 ±0.02. For Af, KcC^v our models predict 
an average value of 2.2btoAi 

X 10^°s ^, and the predicted elastic fraction is 0.77 ± 0.05. 

As noted in I, the work presented in this manuscript can be extended in a number of directions. We can 
apply our model to the description of the semileptonic decays of the light baryons, although these are already 
successfully described by Cabibbo theory. Essentially all experimentally accessible observables for these decays 
have been measured, and it will be interesting to see if our model, constructed with no special reference to chiral 
symmetry or current algebra, can describe the results of these measurements. 

We have not examined the predictions of our model for the many polarization observables which can, in 
principle, be measured in semileptonic decays. In addition, the rare decays of heavy baryons, such as fJfc ^ J7 
can easily be treated in the framework that we have developed. Such processes, along with their meson analogs, 
are used in searches for physics beyond the Standard Model. However, the interpretation of the measured rates 
depend strongly on estimates of the form factors involved (in much the same way that extraction of CKM matrix 
elements depends on the form factors that describe semileptonic decays). Finally, if factorization in some form 
is valid, the semileptonic form factors calculated in the manuscript may also be useful in the description of 
nonleptonic weak decays. 

It may also be possible to systematically improve the quark model used in the present calculation. An obvious 
first step is the implementation of full symmetrization of the spatial wave functions in the Sturmian basis, which 
would allow calculation of results for decays to final state nucleons in this basis. 

We also plan to modify and expand all our baryon spectrum codes to make predictions for baryons containing 
quarks with three different masses. One advantage of this modification is that it will allow us to examine the 
semileptonic decays of Sq. This study will be interesting, as some of the 2n states have an antisymmetric 
(Ag-like) light diquark, while some have a symmetric (fig-like) light diquark p!q |. 
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APPENDIX A: WAVE FUNCTIONS 



The wave function components for states that we consider arc shown here. These components are vahd for 
both the fig and E! states that are treated in this manuscript. For = 1/2+, wave functions are expanded as 



1/2+M 



??f'^V'000000(Pp,PA) +'72'''0OO1OOo(Pp,Pa) +%'^V'000010(Pp,Pa) Xi/2(^^) 



+ ??4"'V'00010i(Pp,Pa)X?/2W [V'lMi,010l(Pp,PA)Xi/2(-^- ^i) 

1/2 



-,1/2 
M 



^"2^^1,0200 

V'2Mz,0002(Pp,PA)xf/2(-^ - ^l) 

For VLq states with = 3/2+, the expansion is 
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(Al) 
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M 

3/2 

M 
3/2 

M 
3/2 

M 

3/2 

M 
3/2 

M 



(A2) 



where [tpLMLnpipnxi\{Pp:Px)xs{M — Ml)Yj^ is a shorthand notation that denotes the Clebsch-Gordan sum 
^j^^{JM\LMl, SM - ML)i)LM^n,i,nxe>,{Pp, Pa)xs(M - Ml). 
For = l/2~ and 3/2~, the expansion is 



ff^Q 111 



V'lMi,0100(Pp, Pa)x5'/2(^ - ^i) 



M 



f2Q 



^i'iMi.oooi(Pp,PA)xf/2(-^ - Ml) 
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M 



(A3) 



where J can take the values 1/2 or 3/2. 
For = 5/2+, the expansion is 



*5/^2+M = "^^fig r?i ''V'2Mx,010l(Pp, Pa)Xv2(^ - Ml) + ■q^'' V'2Mx,020o(Pp, PA)xf/2(^ - Ml) 
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For JP = 



5/2 , the wave function is 




4>nQ ^iMtOooi 



-|5/2 



(A5) 



J-M 



J M 



No other states are expected to have significant overlap with the decaying ground-state Q.q in the spectator 
approximation that we use. 



In this section, we present the analytic expressions we obtained for the form factors, assuming single compo- 
nent wave functions. We list both harmonic oscillator and Sturmian form factors together, ordered by the spin 
and parity of the daughter baryon, beginning with the ground state. We place the form factors for each spin 
and parity in a separate subsection. 

For most of the we treat, there is usually more than one example of the state: the exception is 5/2~, for 
which there is a single state up to A'' = 2. We therefore distinguish among the different states with the same 
by presenting their quark model quantum numbers, in the notation \ns,L,np,ip,n\J.\ >. Here ns takes 
values between 1 and 3, with 1 denoting total quark spin 1/2, with spin wave function of type Xpj 2 denoting 
total quark spin 1/2, with spin wave function of type x\i ^^id 3 denoting total quark spin 3/2. 



APPENDIX B: QUARK MODEL FORM FACTORS 



1. 1/2+, 1200000 > 
a. Harmonic Oscillator Form Factors 
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6. Sturmian Form Factors 
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a. Harmonic Oscillator Form Factors 
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4. 1/2", 1210001 > 
a. Harmonic Oscillator Form Factors 
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